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ABSTRACT 

Context. Extragalactic nuclear activity is best explored with observations at high energies, where the most extreme flux and spectral 
variations are expected to occur, witnessing changes in the accretion flow or in the kinematics of the plasma. In active galactic nuclei 
of blazar type, these variations are the most dramatic. 

Aims. By following blazar outbursts from their onset and by correlating the observed variations at many different wavelengths we can 
reconstruct the behavior of the plasma and map out the development of the flare within the jet. 

Methods. The advent of the Fermi satellite has allowed the start of a systematic and intensive monitoring program of blazars. Blazar 
outbursts are very effectively detected by the LAT instrument in the MeV-GeV domain, and these can be promptly followed up with 
other facilities. Based on a Fermi LAT detection of a high MeV-GeV state, we have observed the blazar PKS 1502+106 with the 
INTEGRAL satellite between 9 and 1 1 August 2008. Simultaneous Swift observations have been also accomplished, as well as optical 
follow-up at the Nordic Optical Telescope. 

Results. The IBIS instrument onboard INTEGRAL detected a source at a position inconsistent with the optical coordinates of 
PKS 1502+106, but consistent with those of the Seyfert 1 galaxy Mkn 841, located at 6.8 arcmin south-west from the blazar, which 
is therefore responsible for all the hard X-ray flux detected by IBIS. At the location of the blazar, IBIS sets an upper limit of ~ 10" 11 
erg s _I cm 2 on the 15-60 keV flux, that turns out to be consistent with a model of inverse Compton scattering accounting for the soft 
X-ray and gamma-ray spectra measured by Swift XRT and Fermi LAT, respectively. The gamma-ray spectrum during the outburst 
indicates substantial variability of the characteristic energy of the inverse Compton component in this blazar. The hard X-ray state 
of the Seyfert appears to be nearly unchanged with respect to the past. On the other hand, its soft X-ray flux (0.3-10 keV) varies 
with a ~50% amplitude on time scales from days to years. On long time scales this is well correlated with the optical flux, with no 
measurable delay. 

Conclusions. In PKS 1502+106 the critical parameters that control variability are the accelerating power transferred to the relativistic 
electrons, and the magnetic field in the emitting region. The spectrum of Mkn 841 in the 0.3-100 keV range is well described by a 
power-law with a cutoff at ~ 150 keV and a Compton reflected continuum. 

Key words, galaxies: active - quasars: individual: PKS 1502+106 - Galaxies: Seyfert - galaxies: individual: Mkn 841 - radiation 
mechanisms: general - gamma-rays: galaxies 



1. Introduction 

The daily survey of the MeV-GeV sky by the Large Area 
Telescope (LAT) instrument (Atwood et al. 2009) onboard the 
Fermi gamma-ray satellite has enhanced dramatically the detec- 
tion rate of flaring blazars, and will allow an accurate determi- 
nation of their variability duty cycle at the highest energies. The 
follow-up of these high states at lower energies and the compar- 
ison with quiescent states make it possible to isolate the basic 
parameters that control the variations at all frequencies and ul- 
timately to set constraints on the behavior of the central engine 
(Ghisellini et al. 2009; Marscher et al. 2010; Abdo et al. 2010a; 
Bottcher & Dermer 2010). 

PKS 1502+106 (z = 1.839, see Abdo et al., 2010b, for a dis- 
cussion on this preferred redshift value) is a bright, variable and 
optically polarized blazar (Impey & Tapia 1988; George et al. 
1994; Lahteenmaki & Valtaoja 1999; Lister & Homan 2005). An 
et al. (2004) suggested, based on the detection of remarkable su- 
perluminal motion with high resolution radio observations, that 
it should be a strong gamma-ray emitter. The source is not in- 
cluded in the Third EGRET Catalog (Hartman et al. 1999), how- 
ever it was detected in the first three months of sky- survey op- 
eration with the LAT (Abdo et al. 2009) and then again with the 
same instrument during a flaring state on 6 August 2008 (Abdo 
et al. 2010b) with a flux above 2 x 10~ 6 photons cirT 2 s~' at 
energies larger than 100 MeV. In January 2009, the source has 
undergone another flare reaching a comparable gamma-ray flux 
level (Abdo et al. 2010b). It is one of the brightest, most intrin- 
sically luminous, and violently variable AGNs of the Fermi LAT 
sample (Abdo et al. 2010c; Abdo et al. 2010d). 

The field of PKS 1502+106 had been observed on 24-26 
January 2006 by INTEGRAL IBIS for a total integration time 
of 83 ks. The observation, analysed by Ubertini et al. (2008) in 
the light of the Fermi report, resulted in the detection of a new 
source, IGR J15039+1022, at a flux level of 1.6 mCrab (i.e., 
2.45 x 10~ n erg cirT 2 s" 1 in the 18-60 keV rangd__). 

The IBIS source, located at RA = 15:03:53.76 and Dec = 
+ 10:22:08.4 (error radius of ~5 arcmin), has been tentatively 
identified with the Seyfert 1 galaxy Mkn 841 (z = 0.036), and 
has been included in the Fourth INTEGRAL IBIS/ISGRI soft 
gamma-ray survey catalog (Bird et al. 2010) and in the second 
INTEGRAL AGN catalog (Beckmann et al. 2009). Ubertini et al. 
(2008) note that PKS 1502+106 is 11 arcmin away from IGR 
J15039+1022 and is not detected in the January 2006 IBIS ob- 
servation. They derive a 2-sigma upper limit flux of 0.7 mCrab 
for the blazar in the range 18-60 keV, assuming a Crab-like spec- 
trum. 

Prompted by the Fermi detection of the gamma-ray out- 
burst, we activated INTEGRAL observations starting on 9 August 
2008 (see preliminary images in Blazek & Hudec 2009). 
Simultaneous observations with Swift XRT and with the Nordic 
Optical Telescope (NOT) were also organized. We report here 
the results and discuss them in the context of AGN models. 



* Based on observations made with: INTEGRAL, an ESA project 
with instruments and Science Data Centre funded by ESA member 
states (especially the PI countries: Denmark, France, Germany, Italy, 
Switzerland, Spain), Czech Republic and Poland, and with the partici- 
pation of Russia and the USA, and the Nordic Optical Telescope, op- 
erated on the island of La Palma jointly by Denmark, Finland, Iceland, 
Norway, and Sweden, in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias. 

1 Note that this flux supersedes the value of 1.2 x 1CT 11 erg cm" 2 s" 1 
reported previously (Ubertini et al. 2008; Abdo et al. 2010b). 



2. Observations and data analysis 

2.1. INTEGRAL 

The source PKS 1502+106 (whose catalogued optical J2000 co- 
ordinates are RA = 15:04:24.98, Dec = +10:29:39.2) was ob- 
served by the INTEGRAL satellite (Winkler et al. 2003) for 203 
ks in Revolution 711 between 2008 August 09, 01:53 UT and 
2008 August 11, 15:12 UT. The on-source time for IBIS/ISGRI 
(Ubertini et al. 2003; Lebrun et al. 2003) was 194.5 ks. The tar- 
get is in the JEM-X (Lund et al. 2003) field of view for 46 ks in 
Rev. 711. 

The screening, reduction, and analysis of the INTEGRAL 
IBIS/ISGRI and JEM-X data have been performed using the 
INTEGRAL Offline Scientific Analysis (OSA) v. 7.0 and v. 8.0, 
respectivelyjjublicly available through the INTEGRAL Science 
Data Centefl (ISDC, Courvoisier et al. 2003). The IBIS/ISGRI 
and JEM-X data have been accumulated into final images 
with algorithms implemented in the software as described in 
Goldwurm et al. (2003) and Westergaard et al. (2003), respec- 
tively. Following similar methods, we have also re-analyzed the 
archival INTEGRAL IBIS and JEM-X observation obtained on 
24-26 January 2006 (Ubertini et al. 2008). 

The source has been also observed with the INTEGRAL 
OMC (Mas-Hesse et al. 2003) with a standard V-band Johnson 
filter. The data have been extracted with default settings, using a 
3x3 pixels binning, which is appropriate for point-like sources. 

2.2. Swift 

X-ray and optical/UV observations of the blazar field have been 
accomplished from 7 to 22 August 2008 with Swift XRT and 
UVOT (Gehrels et al. 2004; Burrows et al. 2005; Roming et 
al. 2005), respectively. The pointing included the location of 
the Seyfert 1 Mkn 841. The results of the analysis of the data 
of PKS 1502+106 have been reported by Abdo et al. (2010b). 
Following similar procedures as those adopted by those authors, 
we have independently reduced the XRT and UVOT data of the 
blazar and of the Seyfert 1. Like in Abdo et al. (2010b), for the 
blazar, given its low count rate (< 0.5 counts s" 1 ), only photon 
counting data were considered and XRT grades 0-12 were se- 
lected. No pile-up correction was required. For the Seyfert, we 
also considered window timing mode data. The source events 
were extracted in circular regions centered on the sources with 
radii depending on the source intensity (10-20 pixels, 1 pixel 
~ 2'.'37), while background events were extracted in source-free 
annular regions around the sources. 

XRT spectra were extracted for each of the two sources 
from the observation of 8 August 2008, that is the one with 
the longest exposure time (12.5 ks). Ancillary response files 
were generated with xrtmkarf and account for different extrac- 
tion regions, vignetting and Point-Spread Function (PSF) cor- 
rections. We used the vOlO spectral redistribution matrices in 
the Calibration Database maintained by HEASARC. All spectra 
were rebinned with a minimum of 20 counts per energy bin to 
allow x 1 fitting within XSPEC (vll.3.2). For the spectral fitting 
we fixed the hydrogen column density to the Galactic value of 
2.22 X 10 20 cm" 2 (Kalberla et al. 2005). 

We extracted XRT light curves from the same image regions 
as for the spectra for source and background in the band 0.2- 
lOkeV. For our analysis, we considered time bins ranging from 
1000 to 12000 s, to ensure that each light curve point would 
have at least 30 source-plus-background counts. The exposure 

2 http://www.isdc.unige.ch/INTEGRAL/download/osa_sw 
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Table 1. Log of Optical Observations at the Nordic Optical 
Telescope 
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* Not corrected for Galactic dust absorption. Errors are 
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requirements were that the bins be at least 50 % exposed. In each 
case, the light curves were corrected for PSF losses, i.e., losses 
due to the extraction region geometry, and bad/hot pixels and 
columns falling within this region. 

Likewise, the UVOT data of both sources were reduced and 
analyzed following procedures similar to those used by Abdo 
et al. (2010b). The photometry of Mkn 841 may include partial 
contamination from the host galaxy bulge. We used an aperture 
of 6" for the v,b,u filter images and 12" for the UV filters. We 
have also reduced and analyzed the Mkn 841 data of 3 February 
2007, that have been included in the Seyfert UVOT study of 
Grupe et al. (2010). 

2.3. Nordic Optical Telescope 

We have observed the IBIS field on 13 August 2008 with the 
ALFOSC camera and optical BVRI filters at the NOT (Canary 
Islands). The data have been reduced following standard pro- 
cedures within the IRAF package (overscan subtraction, bias 
correction, twilight flat fielding). Aperture photometry was per- 
formed with the IRAF routine APPHOT. Eight field stars were 
used to estimate the photometric zero point and their magnitudes 
were taken from the SDSS DR7 (using CDS/Aladin). The field 
stars SDSS gri- magnitudes were converted to BVRI bands fol- 
lowing the transformation^ recommended for the SDSS Data 
Release 5. See Table 1 for a log of the observations. 

3. Results 

3.1. INTEGRAL 

A source is detected in hard X-rays with IBIS/ISGRI at coor- 
dinates RA = 15:03:53.5, Dec = +10:26:45, with an error of 
2.4 arcmin at 90% confidence level (see Figure 1). The detec- 
tion has an 8.6-cr significance in the co-added image of the full 
observation. The centroid of the error circle of this source is 
located at 8.3 arcmin from the optical source corresponding to 
PKS 1502+106, and thus it is inconsistent with the blazar at a 4- 
<t significance level (i.e., the association probability is ~ 10~ 3 ). 

However, it is located at only 1.9 arcmin from Mkn 841 
(optical catalogue coordinates J2000: RA = 15:04:01.20, Dec = 
+ 10:26:16.2), and therefore consistent with it (see Fig. 1). The 
error circle of the IBIS detection of January 2006 is also reported 



3 http://www.sdss.org/DR5/algorithms/sdssUBVRITransform.html 
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Fig. 1. Digitized Sky Survey image of the field of 
PKS 1502+106 and Mkn 841 in red filter. The image ver- 
tical size corresponds to 7 arcmin. The 2 optical sources are 
indicated. Over-plotted are the error circles of the IBIS detec- 
tions of January 2006 (larger circle, with a radius of 5 arcmin) 
and August 2008 (smaller circle, radius of 2.4 arcmin) 



in Fig. 1: as already noted, this hard X-ray source is also consis- 
tent only with the Seyfert galaxy. 

From the count rates measured in our IBIS observation we 
derive a flux of (2.4 + 0.3) xl0~ n erg s _1 crrT 2 in the range 
15-60 keV for Mkn 841, assuming the spectrum measured by 
BeppoSAX (F = 2, Petrucci et al. 2002). The limited signal- 
to-noise ratio does not allow us to extract a spectrum from the 
present IBIS observation), and a 3-cr upper limit of 9.3 x 10~ 12 
erg s^ 1 cm -2 in the same energy range for the blazar. These val- 
ues are similar to the ones determined by Ubertini et al. (2008) 
from the January 2006 IBIS observation, taking into account the 
slight difference in energy range. 

Neither PKS 1502+106 nor Mkn 841 are detected by JEM- 
X in both 2006 and 2008. Assuming a Crab-like spectrum, we 
estimate 3-cr upper limits on the 3-10 keV flux of both sources 
of3.3xl0" n ergs" 1 cirT 2 and 2.3x10-" ergs" 1 cm" 2 for 2006 
and 2008, respectively (consistent with the XRT measurements, 
see below). 

The OMC has measured an optical magnitude of V = 15.4 + 
0.4 for PKS 1502+106. Since this is substantially brighter than 
the V-band flux measured by NOT (see Table 1) and consider- 
ing the relatively large OMC pixel size (about 18 arcsec), we 
have checked for possible contamination from nearby sources. 
Four USNO stars are located within 1 arcmin of the blazar: the 
summed contribution of their magnitudes, when added to the 
flux of the blazar as measured with the NOT, is comparable to 
the measured magnitude so that they contaminate significantly 
the blazar. Therefore, we will ignore its OMC optical measure- 
ment. The OMC field of view included also Mkn 841, for which 
we measured an average magnitude of V = 14.33 ± 0.02 with 
no significant variability. No nearby source contaminates signif- 
icantly the Seyfert optical emission. 

3.2. Fermi LAT 

We have re-analyzed the Fermi LAT spectrum in outburst state 
presented by Abdo et al. (2010b) using 2 different time intervals 
for signal integration: for the same period of five days between 5 
and 10 August 2008 adopted by Abdo et al. (2010b) we obtain a 
photon index of F = 2.05 + 0.04, consistent with that measured 
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by those authors (T = 2.06 + 0.02). For the day of the gamma-ray 
flare peak, 8 August 2008, we find a photon index of T = 1.88 + 
0.08. These indices are reported in Figure 3 (see Discussion). 

3.3. Swift 

Comparison of our results with those of Abdo et al. (2010b) for 
PKS 1502+106 and of Grupe et al. (2010) for Mkn 841 shows 
a general consistency. Swift XRT and UVOT light curves for 
PKS 1502+106 in August 2008 have been reported by Abdo et 
al. (2010b), who found a decreasing trend in X-rays, the after- 
math of the gamma-ray outburst, and a fully resolved flare in 
optical/UV, correlated with the gamma-ray event. In Figure 2a 
we show the light curves of Mkn 841 during the present cam- 
paign and in Figure 2b we report all XRT and UVOT V-band 
flux measurements of the Seyfert galaxy, to compare its long- 
and short-term variability in X-ray and optical. Significant X- 
ray variability of up to ~60% was detected on a time scale of 
10 days and of ~50% in almost 2 years. The V-band variabil- 
ity is less pronounced on both time scales, only about 15-20%. 
The light curves in the other UVOT filters have a behavior sim- 
ilar to the V-band, with somewhat larger variability amplitude, 
although the UV filter data have a larger error than the optical 
ones. It is interesting to note that the long-term X-ray and opti- 
cal variations are correlated (Fig. 2b). 

The XRT spectrum of the blazar on 8 August 2008 is satis- 
factorily (x 2 = 32.5 for 28 d.o.f.) fitted by a single power-law 
Fe °c E~ r with photon index r = 1 .52 + 0.09 (errors on fitted 
X-ray spectral parameters are at 90% confidence), and the flux in 
2-10keV is 1.8x 10~ 12 erg s _1 cirT 2 . These results are consistent 
with those of Abdo et al. (2010b), who integrated the X-ray spec- 
tral signal over 4 days around the outburst. The extrapolation of 
the XRT flux to the IBIS energy range of 15-60 keV yields a flux 
of 3.8 x 10~ 12 erg s _1 cirT 2 , which is 6 times lower than the flux 
measured by IBIS for the detected source. This independently 
suggests that, barring an improbably large spectral change be- 
tween the soft and hard X-ray ranges, the IBIS source does not 
coincide with the blazar. 

A single power-law describes well (x 2 = 46.3 for 52 d.o.f.) 
the spectrum of Mkn 841 of 8 August 2008. The photon index 
is F = 2.03 + 0.07 and the fitted flux in 2-10 keV is 1.2 x 10" 11 
erg s _1 crrT 2 . The extrapolation to the IBIS energy range yields 
a flux of 9.9 x 10~ 12 erg s _1 cnr 2 (15-60 keV), which is 2.4 
times lower than measured by IBIS, and therefore compatible 
with it, both taking into account the not exact simultaneity of 
the XRT spectrum and IBIS detection and assuming a spectral 
change between soft and hard X-ray range due to the emergence 
of a reflected Compton hump (see Discussion). No significant 
X-ray spectral variability is seen during the August 2008 cam- 
paign. Comparing our XRT spectrum of Mkn 841 with those of 
previous epochs reported by Grupe et al. (2010) we find some 
modest, marginally significant X-ray spectral variability, corre- 
lated with the flux in the sense of a harder spectrum in lower 
X-ray emission state. The broad-band UVOT spectra indicate a 
rather constant UV-optical spectral slope. 

3.4. Nordic Optical Telescope 

The results of the NOT observations are reported in Table 1. 
For the photometry of Mkn 841 we used a rather small aper- 
ture (5 pixel radius equivalent to 0.95 arcsec) in order to reduce 
the host galaxy contribution. However, we noted that the profile 
of one of the comparison star and the target are very similar in B- 



and V-band images, leading us to conclude that the host galaxy 
contribution is negligible. The NOT measurement in V-band of 
Mkn 841 is consistent with the INTEGRAL OMC measurement 
within the errors and taking into account the difference in filter 
central wavelength. 

4. Discussion 

4.1. PKS 1502+106 

The data of our multiwavelength campaign in August 2008 are 
reported in Figure 3, together with archival non-simultaneous 
data that sample a quiescent broad-band spectral state. Also 
shown is the Fermi LAT spectrum detected during the first 
3 months of the survey (Abdo et al. 2009). Abdo et al. 
(2010b) have modeled the SED of PKS 1502+106 with a syn- 
chrotron component at the radio-to-UV wavelengths and inverse 
Compton scattering at the higher energies. In particular, syn- 
chrotron self-Compton radiation is mostly responsible for the X- 
ray spectrum, while external Compton upscattering is invoked to 
reproduce the MeV-GeV spectrum that has a curved shape in the 
lower state and a power-law slope during the flare (5-10 August 
2010). 

In order to fit both the August 2008 low and high states 
and the archival SEDs we applied a simple, one-zone, ho- 
mogeneous leptonic model based on synchrotron radiation at 
radio-to-UV frequencies, and inverse Compton upscattering of 
both synchrotron photons (self-Compton) and external photons 
of the broad emission line region (BLR), dominating at X-ray 
and gamma-ray frequencies, respectively, following an approach 
similar to that of Abdo et al. (2010b). The details of the model 
are described in Ghisellini, Tavecchio, & Ghirlanda (2009) and 
Ghisellini et al. (2010). A strong constraint of the model is 
the variation of the MeV-GeV spectrum in August 2008, indi- 
cating that the peak frequency of the external Compton com- 
ponent changes significantly from low to outburst state. Note 
that the INTEGRAL IBIS upper limit also sets a stringent con- 
straint. Although the optical-UV spectrum follows a steep single 
power-law and does not require an additional thermal compo- 
nent (Fig. 3), the presence of an external Compton component at 
gamma-rays demands a powering source of the BLR. We iden- 
tify this with an accretion disk whose luminosity is limited by 
the smooth shape of the optical-UV spectrum to be no larger 
than L d ^ 1.7 x 10 46 erg s _1 , consistent with that postulated by 
Abdo et al. (2010b). This is ~11% of the Eddington luminosity 
for an assumed central black hole mass of ~ 1O 9 M , therefore 
in this blazar the jet power exceeds by a large amount the accre- 
tion power (see Ghisellini et al. 2010). The radius of the BLR is 
assumed to be Rblr = 4.1 X 10 17 cm. The total X-ray corona 
luminosity is assumed to be 30% of L&, its spectral shape is as- 
sumed to be oc v _1 exp(-/zv/150 keV). The model parameters for 
the present SEDs are reported in Table 2 and the corresponding 
model curves are overplotted to the data in Figure 3 in the same 
colors indicated in the table. The jet viewing angle 8 V is 2.3° for 
all states. 

Our SED analysis differs from that of Abdo et al. (2010b) 
in the changes of the parameters that are responsible for vari- 
ability among the various emission states (see Table 2). At vari- 
ance with Abdo et al. (2010b) the parameters that exhibit the 
largest variations in our modeling are: the magnetic field B, that 
increases by an order of magnitude between the highest and low- 
est gamma-ray state; the total power injected in the jet P'., that 
is a factor 10 higher in the highest than in the lowest gamma-ray 
state; the random Lorentz factor jb of the electron distribution 
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Fig. 2. Swift light curves of Mkn 841 : (a) in August 2008 and (b) in 2007-2008: XRT count rates (in c s" 1 , open circles) and UVOT 
fluxes (filled circles) obtained with an arbitrary magnitude-to-flux conversion coefficient from the V-band magnitudes. The errors 
are statistical only and correspond to 1-cr. The data points in (b) are from Grupe et al. (2010), except the last XRT and UVOT ones, 
that correspond to the averages of the August 2008 points for the two instruments shown in (a). The zero time in (b) corresponds to 
2007 January 1,00 UT 



cooling break, that is significantly larger in the highest gamma- 
ray state, consistent with the harder spectrum. The variation of 
the magnetic field is consistent with the dominance of the in- 
verse Compton cooling in the brightest gamma-ray state of 8 
August 2008: a higher intensity of the inverse Compton scatter- 
ing radiation implies that the ratio between the radiation energy 
density and the magnetic field energy density is relatively large, 
while the magnetic energy density dominates this ratio in lower 
inverse Compton states. The increase of the injected power P'- 
is related to an enhancement of the activity of the acceleration 
mechanism, so that the observed correlation between P' i and y\, 
is not unexpected. We also find that the shape of the particle 
distribution changes between the low and the hard states, with 
harder slopes below jb and steeper above. The synchrotron com- 
ponent, traced by the optical data, does not vary as dramatically 
as the inverse Compton component (Me V-GeV) in the same time 
scale, a behavior seen, for instance, in the blazars 3C 279 (see 
Wehrle et al. 1998 and Giuliani et al. 2009) and PKS 1454-354 
(Ghisellini et al. 2009), but different from that of other well 
monitored blazar sources like 3C 454.3 (Bonning et al. 2009; 
Vercellone et al. 2009; Donnarumma et al. 2009; Vercellone et 
al. 2010; Pacciani et al. 2010) and PKS 1510-089 (Abdo et al. 
2010e; D' Ammando et al. 2010), where the correlated variations 
in optical and gamma-rays have comparable amplitudes on time 
scales from days to years (see however Ghisellini et al. 2007). 

A multiwavelength monitoring of few bright Fermi blazar 
sources on a range of time scales from days to few years would 
allow us to obtain sensitive insight into the undoubtedly com- 
plex jet behavior. If used systematically in these campaigns, 
INTEGRAL will sample, or will set constraining limits on, the 
flux at hard X-ray energies, that is crucial for locating precisely 
the peak of the high energy component. Besides this intrinsic 
merit, INTEGRAL IBIS - with its large field-of-view - often al- 
lows serendipitous detections of sources located within a few 
degrees of the primary target. Both in the present case, with the 
detection of the Seyfert Mkn 841, and in other previous circum- 



stances (e.g., Pian et al. 2005), this advantage has led to an ad- 
ditional contribution in the study of the hard X-ray extragalactic 
sources. 

4.2. Mkn 841 

BeppoSAX observations of the Seyfert 1 Mkn 841 in 1999 and si- 
multaneous XMM-Newton and BeppoSAX observations in 2001 
showed a high energy cut-off at around 100 keV (Bianchi et al. 
2001; Petrucci et al. 2002; Bianchi et al. 2004). As noted by 
Ubertini et al. (2008), this strongly suggests that this source, 
despite being 20 times brighter in X-rays than PKS 1502+106 
when quiescent, is unlikely to emit in the Fermi LAT domain at 
the reported level (Abdo et al. 2010b). Therefore, the Fermi LAT 
detection can be firmly identified with emission from the blazar. 
On the other hand, the imaging and spectral arguments made 
in the previous section suggest that the hard X-rays detected by 
IBIS should be entirely attributed to the Seyfert. 

From the January 2006 and present INTEGRAL IBIS ob- 
servations there is no evidence of significant variability of the 
hard X-ray flux of Mkn 841. This source was first detected in 
hard X-rays (50-150 keV) in 1996 by CGRO OSSE (Zdziarski, 
Poutanen, & Johnson 2000), at a flux level of (6 ± 2) x 10" 11 
erg s _1 cm 2 , an order of magnitude higher than the extrapola- 
tion of the IBIS flux in the OSSE band using a reflected Compton 
spectrum (see below and Fig. 4), a difference that is not to be 
considered dramatic owing to the low significance of the OSSE 
detection (we tend to exclude contamination of the OSSE detec- 
tion by PKS 1502+106, see below). On the other hand, the Swift 
BAT detection of Mkn 841 obtained during the first 39 months of 
the mission (Cusumano et al. 2010) is somewhat lower than the 
IBIS flux (see Figure 4). This indicates a moderate secular vari- 
ability of the hard X-ray emission. The softer (0.3-10 keV) X-ray 
flux varies more substantially during the Swift mission lifetime, 
and exhibits maximum amplitude variability of about ~50% on 
both months and days time scales (Fig. 2), in agreement with the 
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behavior observed by other satellites, that detected even faster 
variations (George et al. 1993; Nandra et al. 1995; Bianchi et al. 
2001;Petrucci etal. 2007). 

In Figure 4 we report the spectral energy distribution of 
Mkn 841 based on the INTEGRAL, Swift, and NOT quasi- 
simultaneous (within few days) observations of August 2008, 
and compare it with the historical data and with the Swift XRT 
and UVOT observations of February 2007. The Swift BAT aver- 
age flux and spectrum from the 39-month survey (Cusumano et 
al. 2010) are also shown. The spectrum has the overall typical 
shape of a Seyfert 1 nucleus, with a prominent rising optical/UV 
component, a relatively soft spectrum in the 0.3-10 keV range 
and a hard X-ray component with a cutoff around 150 keV, that 
suggests a relatively high coronal maximum temperature. The 
NOT, UVOT and XRT measurements of August 2008 indicate a 
relatively high state. The limited accuracy of XRT does not allow 
us to appreciate the possible presence of the soft X-ray excess, 
that is often observed in Seyferts and has been first detected pre- 
cisely in this object (Arnaud et al. 1985). However, as already 
noted, the lower state of February 2007 is also harder, suggest- 
ing a variable soft component, more prominent in the high states. 
The correlation between X-ray spectral index and flux (harder X- 
ray spectrum during lower X-ray states) is confirmed by all Swift 
XRT observations of this source (see Grupe et al. 2010). 

We have overimposed to the data, without any scaling, the 
model of a black body with characteristic temperature of ~100 
e V to reproduce phenomenologically the soft X-ray excess, plus 
a power-law due to Compton upscattering of the coronal soft X- 
rays, with a cutoff at 100 keV, and a reflected continuum and 
fluorescent iron line component at the hard X-rays, proposed 
by Bianchi et al. (2001) as the best description of the 1999 
BeppoSAX spectrum (see also Ponti et al. 2008 for a Suzaku 
observation of the hard X-ray spectrum). Apart from the iron 
line, that is out of the reach of the limited sensitivity and spec- 
tral resolution of XRT, the model reproduces satisfactorily the 
XRT, IBIS and BAT data (Fig. 4). 

The average optical-to-X-ray colour index a ox obtained from 
the present and previous (Grupe et al. 2010) XRT and UVOT 
observations of Mkn 841 is consistent with the relationship be- 
tween a ox and UV luminosity found in optically selected AGNs 
(Strateva et al. 2005; Steffen et al. 2006), although this is not 
the case when the a ox values of the individual Swift observations 
are compared, because the X-ray variations are larger than the 
optical ones. Correlated X-ray and optical variability with larger 
amplitude in X-rays than optical has been seen in other Seyferts 
(e.g., Breedt et al. 2010), although the opposite is also seen in ob- 
jects of this AGN class. This behavior of Mkn 841 optical-X-ray 
variability and the correlated X-ray flux and spectral variations 
(steeper spectrum for brighter states) must be both taken into ac- 
count when constructing an emission scenario for this source. 
However, our data are not adequate to provide a detailed physi- 
cal characterization of the soft X-ray excess component and its 
relation to the optical-UV bump. This requires a complex model 
(see e.g. Petrucci et al. 2007; Ponti et al. 2008; Longinotti et al. 
2010) and an accurate and sensitive simultaneous monitoring in 
optical/UV and X-rays. 

We conclude with a note on the issue of possible contam- 
ination of the hard X-ray flux of Mkn 841 in the observations 
that preceded the INTEGRAL one: the model we used to repro- 
duce the broad-band spectrum of PKS 1502+106 (Fig. 3) pre- 
dicts an intensity of about 5 x 10" 12 erg s cirT 2 in hard X-rays, 
with little variability among very different MeV-GeV emission 
states. This is a rather negligible fraction (20%) of the IBIS de- 
tected flux. Assuming therefore that neither the blazar nor the 



Seyfert vary dramatically in hard X-rays, both the BeppoSAX 
PDS measurement of Mkn 841 (Bianchi et al. 2001; Petrucci et 
al. 2002), which is consistent with IBIS, and the CGRO OSSE 
measurement (Johnson et al. 2000) do not contain significant 
contribution from the blazar, which lies in the field of view of 
the two instruments during Mkn 841 observations. The uncer- 
tainty circle of the Swift BAT Mkn 841 detection has a radius of 
2.5 arcmin (Cusumano et al. 2010), which excludes contamina- 
tion from PKS 1502+106. 
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Fig. 3. Spectral energy distributions of PKS 1502+106. The red circles and the darkgreen bowtie at GeV energies represent the data 
of the outburst state of 8 August 2008; the red bowtie represents the Fermi LAT outburst spectrum averaged over the days 5-10 
August and the downward pointing arrow is the INTEGRAL IBIS/ISGRI upper limit of 9-11 August. The blue circles represent 
the post-flare state of 15-22 August 2008 (Swift XRT and UVOT) and the blue bowtie is the spectrum in the low gamma-ray state 
observed by Fermi during its first 3 months in orbit (Abdo et al. 2009), equivalent to the Fermi LAT post-flare state after 1 1 August 
2008 (see Abdo et al. 2010b). The NOT data (red open squares) refer to 13 August 2008. The ASCA and CGRO EGRET data (blue) 
and the grey open circles represent a historical quiescent radio-to-gamma-ray state. The optical/UV data have been corrected for 
dust extinction in our Galaxy using Eb-v = 0.032 (Schlegel et al. 1998) and the extinction curve of Cardelli et al. (1989). The 
optical magnitudes from NOT and the optical-UV magnitudes from UVOT have been converted to fluxes according to Fukugita, 
Shimasaku, & Ichikawa (1995) and to Poole et al. (2008), respectively. The errors include the statistical and systematic uncertainties. 
The solid curves represent the synchrotron plus synchrotron self-Compton and external Compton model we used to fit the data of the 
highest gamma-ray state (light green), average outburst gamma-ray state (orange), the low state (cyan) and the historical quiescent 
state (black). See text and Table 2 for the model parameters. The black dashed curve represents the thermal emission from the dust 
torus, the accretion disk and its corona. The models do not account for the radio emission, that is likely produced in larger regions 
than the homogeneous one where the high energy emission takes place. 
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Fig. 4. Spectral energy distributions of Mkn 841 in August 2008 (red symbols: INTEGRAL IBIS/ISGRI, 9-11 August; Swift XRT 
spectrum, 8 August; Swift UVOT, 13-15 August; NOT photometry, 13 August), on 3 February 2007 (green symbols: Swift XRT and 
UVOT) and at prior epochs, not from simultaneous data (taken from the NED Extragalactic Database, grey squares). The optical/UV 
data have been corrected for dust extinction in our Galaxy using Eb-v - 0.03 (Schlegel et al. 1998) and the extinction curve of 
Cardelli et al. (1989). The optical magnitudes from NOT and the optical-UV magnitudes from UVOT have been converted to fluxes 
according to Fukugita, Shimasaku, & Ichikawa (1995) and to Poole et al. (2008), respectively. In addition, the Swift BAT spectrum 
from the first 39 months of the mission is shown in blue (December 2004 - February 2008, Cusumano et al. 2010). The errors include 
the statistical and systematic uncertainties. The model fitted by Bianchi et al. (2001) to the 1999 BeppoSAX spectrum, including 
a soft X-ray excess, a Comptonized spectrum of a hot corona with a high energy cutoff and a reflected component, is shown as a 
dashed curve. 
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Table 2. Parameters of the theoretical SEDs of PKS 1502+106. 
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* dissipation radius in units of 10 15 cm and (in parenthesis) in units of gravitational radii 
c power injected in the blob calculated in the comoving frame, in units of 10 45 erg s _I 
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